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Abstract. We develop a linear systems theory that coincides with the ex- 
isting theories for continuous and discrete dynamical systems, but that also 
extends to linear systems defined on nonuniform time scales. The approach 
here is based on generalized Laplace transform methods (e.g. shifts and con- 
volution) from the recent work |13) . We study controllability in terms of the 
controllability Gramian and various rank conditions (including Kalman's) in 
both the time invariant and time varying settings and compare the results. 
We explore observability in terms of both Gramian and rank conditions and 
establish related realizability results. We conclude by applying this systems 
theory to connect exponential and BIBO stability problems in this general 
setting. Numerous examples are included to show the utility of these results. 



1. Introduction 

In this paper, our goal is to develop the foundation for a comprehensive linear 
systems theory which not only coincides with the existing canonical systems theories 
in the continuous and discrete cases, but also to extend those theories to dynamical 
systems with nonuniform domains (e.g. the quantum time scale used in quantum 
calculus [9]). We quickly see that the standard arguments on R and Z do not go 
through when the graininess of the underlying time scale is not uniform, but we 
overcome this obstacle by taking an approach rooted in recent generalized Laplace 
transform methods [131 119j . For those not familiar with the rapidly expanding area 
of dynamic equations on time scales, excellent references are [SJ [7] . 

We examine the foundational notions of controllability, observability, realiz- 
ability, and stability commonly dealt with in linear systems and control theory 
[31 [221 El] • Our focus here is how to generalize these concepts to the nonuni- 
form domain setting while at the same time preserving and unifying the well-known 
bodies of knowledge on these subjects in the continuous and discrete cases. This 
generalized framework has already shown promising application to adaptive control 
regimes [TB I fT7 ] . 
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Throughout this work, we assume the following: 

• T is a time scale that is unbounded above but with bounded graininess, 

• A(t) € R nxn ,B(t) € R nxm ,C(t) £ RP xn , and D(t) G RP xm are all rd- 
continuous on T, 

• all systems in question are regressive. 

The third assumption implies that the matrix / + n(t)A(t) is invertible, and so 
on Z, the transition matrix will always be invertible. We are therefore justified in 
talking about controllability rather than reachability which is common [SJ Q31 [25] 
since the transition matrix in general need not be invertible for T = Z. 

In the following sections, we begin with the time varying case, and then proceed 
to treat the time invariant case. We will get stronger (necessary and sufficient) 
results in the more restrictive time invariant setting relative to the time varying 
case (sufficient conditions) . Although some of the statements contained in this work 
can be found elsewhere [U [15] , in each of these cases proofs are either missing, 
are restricted to time invariant systems, or are believed to be in error [15] when 
T has nonuniform graininess. Moreover — and very importantly — the methods used 
here are rooted in Laplace transform techniques (shifts and convolution) , and thus 
are fundamentally different than the approaches taken elsewhere in the literature. 
This tack overcomes the subtle problems that arise in the arguments found in [15] 
when the graininess is nonconstant. 

2. Controllability 

2.1. Time Varying Case. In linear systems theory, we say that a system is con- 
trollable provided the solution of the relevant dynamical system (discrete, continu- 
ous, or hybrid) can be driven to a specified final state in finite time. We make this 
precise now. 

Definition 2.1. Let A{t) e R" x, \ B(t) £ R nxm , C(t) e W xn , and D(t) E W xm 
all be rd-continuous matrix functions defined on T, with p,m < n. The regressive 
linear system 

x A (t) = A(t)x{t) + B(t)u(t), x(t ) = x , 

y(t) = C(t)x(t)+D(t)u(t), 

is controllable on [to,tf] if given any initial state xo there exists a rd-continuous 
input u(t) such that the corresponding solution of the system satisfies x(tf) = Xf. 

Our first result establishes that a necessary and sufficient condition for con- 
trollability of the linear system (|2.ip is the invertibility of an associated Gramian 
matrix. 

Theorem 2.2 (Controllability Gramian Condition). The regressive linear system 

x A (t) = A(t)x{t) + B(t)u(t), x(t Q ) = x , 

y(t) = C{t)x(t)+D(t)u(t), 

is controllable on [io,i/] if and only if the n x n controllability Gramian matrix 
given by 

Sfc(to, t f ) := I ' ^ A (t ,(r(t))B(t)B T (t)^l(t , a(t)) At, 

is invertible, where <&z{t,to) is the transition matrix for the system X A (t) = 
Z(t)X(t),X(t )=L 
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Proof. Suppose ^c(to,t/) is invertible. Then, given xq and Xf, we can choose the 
input signal u(t) as 

u(t) = -B T (t)^(t ,a(t))^ 1 (t ,t f )(x - <f> A (t ,tf)x f ), f G [to,*/), 

and extend u(t) continuously for all other values of i. The corresponding solution 
of the system at * = */ can be written as 

x(t f ) = $ A (t f ,t )x + f * $ A (t f ,a(t))B(t)u(t)At 

= ®A(tf,t )x 

- / $ A {tf,a(t))B(t)B T (t)$l(t ,a(t))&c\to,tf)(x -<!>A{t ,t f )x f )At 

Jto 

= $ A (tf,t )x - $a(*/,*q) 

x f ' $ A {h,a(t))B(t)B T (t)$l(t a ,a(t)) At&^itoit^ixo ~ $A(to,t f )x f ) 

Jt 

= ®A{tf,to)XQ ~ (®A(tf,t )x ~ X f ) 
= x f> 

so that the state equation is controllable on [to,*/]. 

For the converse, suppose that the state equation is controllable, but for the sake 
of a contradiction, assume that the matrix ^c(*o,*/) is not invertible. If £^c(to,*/) 
is not invertible, then there exists a vector x a ^ such that 

= x^ c (to,tf)x a = f 1 xZ$ A (t ,a(t))B(t)B T (t)$l(t , a{t))x a At 

Jt 

= r \\x T a <S> A (t ,a(t))B{t)\\ 2 At, (2.2) 

Jto 

and hence 

xl$> A {t ,a{t))B{t) = Q, *G[* ,*/). (2.3) 

However, the state equation is controllable on [to,t/], and so choosing xq — x a + 
<fr A (tQ,tf)xf, there exists an input signal u a {t) such that 

x f = $ A (t/, t )a; + f 1 $ A (*/, a{t))B{t)u a {t) At, 

Jto 

which is equivalent to the equation 

Xa 

/ <Z> A (t Q ,o-{t))B(t)u a (t)At. 

Jt 

Multiplying through by x^ and using (|2.2p and (|2.3[) yields x^x a = 0, a contradic- 
tion. Thus, the matrix ^c(*o,t/) is invertible. □ 

Since the controllability Gramian is symmetric and positive semidefinite, Theo- 
rem [52] can be interpreted as saying (|2.1[) is controllable on [to,*/] if and only if 
the Gramian is positive definite. A system that is not controllable on [*o,t/] may 
become so when either */ is increased or to is decreased. Likewise, a system that 
is controllable on [to,*/] may become uncontrollable if to is increased and/or t/ is 
decreased. 
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Although the preceding theorem is strong in theory, in practice it is quite lim- 
ited since computing the controllability Gramian requires explicit knowledge of the 
transition matrix, but the transition matrix for time varying problems is generally 
not known and can be difficult to approximate in some cases. This observation 
motivates the following definition and our next theorem. 

Definition 2.3. If T is a time scale such that fj, is sufficiently differentiable with the 
indicated rd-continuous derivatives, define the sequence of n x m matrix functions 

K (t) := B(t), 

K j+1 (t) := (I + rtaVfiAWt)))- 1 *?® [V + mM*))^*)))- 1 (^ A (t)A(a(t)) 

+ Kt)A A (t))(I + rit)A(t))- 1 + A(t)(I + ^(t))" 1 ] 
j=0,l,2,.... 

A straightforward induction proof shows that for all t,s, we have 
d j 

^- l<P A (a(t),o-(s))B( S )} = $ A (cr(t), a(a))Kj{a), j = 0, 1,2, . . . 

Evaluation at s = t yields a relationship between these matrices and those in 
Definition [ 



^(t) = ^j[^A(o-(t),a(s))B(s)] 



J = 0,1, 2,. 



s=t 

This in turn leads to the following sufficient condition for controllability. 

Theorem 2.4 (Controllability Rank Theorem). Suppose q € Z + such that, for 
t € [to,tf], B(t) is q-times rd- continuously differentiable and both of fi(t) and A(t) 
are (q — l)-times rd- continuously differentiable. Then the regressive linear system 

x A (t) = A(t)x(t) + B(t)u(t), x(t Q ) = x Q , 

y(t) = C(t)x(t) + D{t)u(t), 

is controllable on [to,tt] if for some t c € [to,tt), we have 

mnk[K {t c ) Ki(t c ) ... K q (t c )] = n, 

where 

AsJ 



K J (t) = — [^ A (a(t),a( S ))B( S )} 



J = 0, 1, . 



Proof. Suppose there is some t c £ [to,tf) such that the rank condition holds. For 
the sake of a contradiction, suppose that the state equation is not controllable on 
[to,tf). Then the controllability Gramian &c(to, tf) is not invertible and, as in the 
proof of Theorem 12.21 there exists a nonzero n x 1 vector x a such that 

x T a § A {ta, o-{t))B(t) =0, t G [to, t f ). (2.4) 

If we choose the nonzero vector Xb so that Xb = ^(io, o-(t c ))x a , then (12. 4p yields 

x£&A((r(tc), o{t))B{t) =0, t G [to, i/). 

In particular, at t — t c , we have x^ Ko(t c ) = 0. Differentiating ()2.4|) with respect 
to t, 

xl^ A {o-{t c ),a[t))Kx{t) = 0, t G [to,*/), 
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so that xjK\{t c ) = 0. In general, 
d? 



d? 

— [ x Z<S>l(a(t c ),a(t))B(t)] ^=xlK j {t c ) = Q, j = 0, 1, 



Thus, 

xl[K (t B ) Jf x (t c ) ... ^(i c )]=0, 

which contradicts the linear independence of the rows guaranteed by the rank con- 
dition. Hence, the equation is controllable on [£o,t/]- D 

When T = M, the collection of matrices Kj (t) above is such that each member is 
the jth derivative of the matrix ^A(o'(t),a(s))B(s) = $a(£> s)B(s). This coincides 
with the literature in the continuous case (see, for example, [3[8l[24]). However, 
while still tractable, in general the collection Kj(t) is nontrivial to compute. The 
mechanics are more involved even on Z, which is still a very "tame" time scale. 
Therefore, the complications of extending the usual theory to the general time 
scales case are evident even at this early juncture. 

Furthermore, the preceding theorem shows that if the rank condition holds for 
some q and some t c € [to,tf), then the linear state equation is controllable on 
any interval [io,f/] containing t c . This strong conclusion partly explains why the 
condition is only a sufficient one. 

2.2. Time Invariant Case. We now turn our attention to establishing results 
concerning the controllability of regressive linear time invariant systems. The gen- 
eralized Laplace transform presented in [13j [JJ] allows us to attack the problem in 
ways that simply are not available in the time varying case. 

First we recall a result from DaCunha in order to establish a preliminary technical 
lemma. 

Theorem 2.5. [12] For the system X A (t) = AX(t), X(t Q ) = I, there exist scalar 
functions {^o(t,to), . . . ,7 n _i(i,£o)} C C£?(T,K.) such that the unique solution has 
representation 

n-l 

e A {t,t ) = J2 ^S(Mo)- 
Lemma 2.6. Let A, B E R nx ™ and u :=u Xo (t f , a(s)) e C rd (T, M" xl ). Then 

span | £ f e A (s,t )Bu XQ (t f ,a(s))As} = span{£?, AB, . . . , A'^B}. (2.5) 



1 1 

1 n-l 



Proof. Let {7fc(Mo)}fc=o ^ e ^ ne collection of functions that decompose the expo- 
nential matrix as guaranteed by Theorem 12.51 This collection forms a linearly 
independent set since it can be taken as the solution set of an n-th order sys- 
tem of linear dynamic equations. Apply the Gram-Schmidt process to generate an 
orthonormal collection {jk(t, to)}^Zo- The ^ w0 collections are related by 

[7o (Mo) 7i (Mo) ••• 7,1-1 (Mo)] 



[7b(Mo) 7i(Mo) ••• 7n-i(Mo)] 



Pll Pl2 ■■■ Pin 
P22 ■■ ■ P2n 
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where the matrix on the right is the triangular matrix obtained from the QR fac- 
torization of the vector consisting of the functions {^k{t, io)}fc=o 011 ^ ne ^ e ^- 
Using the QR factorization, we can write the matrix exponential as 



e A {t,t ) =J2-f k {t,t )A h 



k=0 
n-1 



7n-l(Mo)] PkA k , 



k=0 



where pk is the fc-th column vector of the triangular matrix R. It is worth recalling 
here that the entries on the diagonal of this matrix are norms of nonzero vectors 
and are thus positive. That is, pa > for all i. 
Rewriting the integral from (|2.5[) , 

rt f 

e A (s 7 t )Bu Xo (t f ,a(s))As 

rtf 



J2lk{s,t )A k Bu Xo {tf,a{ S ))As 
-J fc=0 
n-l „ tf 

= Y, AkB j k (s,t )u Xo (t f ,<T(s))As 

k=0 Jt ° 

= Y. AhB / [7o(«.*o) 7i(Mo) ••■ 



Let 



Vk 



%-l(s,t )] pkU Xo (tf,a(s)) As. 



[7o(S)*o) 7i(Mo) ••• 7n-i(Mo)] PkU Xo (tf,a(s)) As 



fc = 0, 1, . . . , n — 1. We will show that span{y , Ui, ■ ■ ■ , Un-i} = K n . That is, there 
exists some w S C r d(T, R™ x 1 ) so that for any arbitrary but fixed collection of vectors 
{z 0l Z\,..., z„-i} C IR™* 1 , the system 



Pll%(s,t )u Xo (tf,a(s))As : = z = 



la 



zoo 
zoi 

z 0(n~l). 



(7o(Mo)Pi2 +li{s,t )p 2 2)u Xo {tf,(j{s))As := z\ 



Z\0 



z \{n-\) 



(7o(Mo)Pln H r-7n-l(S)*o)P (t f ,a(s))As 

• — z n — i 



Z{n-l)0 
z (n-l)l 

z (n-l)(n-l) 



EJDE-2009/37 



CONTROLLABILITY, OBSERVABILITY, REALIZABILITY 



7 



has a solution. 

To accomplish this, we use the fact that the collection jk{s,to) is orthonormal 
and search for a solution of the form 

Ux (tf><r( s )) = ( u j) = /%ii( s > to^j ■ 
Starting with uo, the equations become 

zoo 



r 'If ~ 

/ ToPn As = 

Jt o i=0 

/ (70P12 + 71P22) } y $ii As = zi 

J to „-_n 



/tf ™ _1 
{loPln + 7lP2n H h 7n-lPnn) X! As = z ("-l)0- 

a i=0 

Since the system 7^ is orthonormal, we can simplify the equations above using the 
fact that the integral of cross terms Ji'fj, i 7^ j, is zero. After doing so, the system 
becomes a lower triangular system that can be solved by forward substitution. (The 
observation that pa ^ is crucial here, since this is exactly what allows us to solve 
the system.) For example, the first equation becomes 

/ 7oPn 22 /3°7i As = / to 2 /?SpiiAs = 8°p n = z 01 , 

Jto i=0 Jto 

so that 8® = Using this value for /?§ in the second equation, 



t i = Q 



(70P12 + 71P22) tfili As = / (70P12 + 71P22) {Polo + /3?7i) A 



s 



Pl2 1 «o 

= Z i+ l\V22 

Pll 

— z 10, 

so that 8? = —Zu El2 —zoi- We can continue solving the system in like manner 

^ L P22 P11P22 J 

by using forward substitutions to find 8® for all j — 0, 1, . . . ,n — 1, which will in 

turn yield uo — Y^i=a Ptli- Repeating this process for ux, U2, ■ ■ ■ , u n -i, we find the 
correct linear combinations of -jk to solve the system, and so the claim follows. □ 

We are now in a position to establish the following analogue of the Controllability 
Rank Theorem (Theorem 12. 4p. 

Theorem 2.7 (Kalman Controllability Rank Condition). The time invariant re- 
gressive linear system 

x A (t) = Ax(t) + Bu(t), x(t )=x , 
y(t) = Cx{t)+Du(t), 
is controllable on [to,i/] if and only if the n x nm controllability matrix 

[B AB ■■■ A n-1 5] 



s 



J. M. DAVIS, I. A. GRAVAGNE, B. J. JACKSON, R. J. MARKS II 



E.7DE-2009/37 



has rank n. 

Proof. Suppose the system is controllable, but for the sake of a contradiction that 
the rank condition fails. Then there exists annxl vector x a such that 

x%A k B = 0, ft = 0,...,n-l. 

Now, there are two cases to consider: either x^Xf — or x^Xf ^ 0. 
Suppose x^Xf ^ 0. Then for any t, the solution at time t is given by 

x(t)= / e A (t, a(s))Bu Xo (s) As + e A (t,t )x 
Jt 

= e A (t, 0) * Bu{t) + e A (t, 0)x 

= Bu(t)*e A (t,0) + e A {t, 0)x 
t 

e A (s, t )Bu X0 (t, a (a)) As + e A (t, t )x , 



where we have written the solution as a (time scale) convolution and appealed to 
the commutativity of the convolution [13l[T9j. Choose initial state xo = By, where 
y is arbitrary. Then, again by commutativity of the convolution and Theorem [231 



x a x(t)=x a / e A (s,t )Bu Xo (t,a(s)) As + x a e A (t,to)xo 
J to 

= / ^ k ( s ' A k Bu X0 (t 1 a{s)) As + ^ 7fc (i, t )^A fc Sy 
k=a k=a 

= 0, 

so that x^x(t) = for all t. This is a contradiction since x^x{tj) = x^x/ ^ 0. 

Now suppose x^Xf = 0. This time, we choose initial state xq = (i/, to)x a . 
Similar to the equation above, 

/{ n-i 
^] Jk(s, t Q )xlA k Bu Xo (t, a(s)) As + x^e A {t, t Q )e A ^(t f , t )x a 

= xJe A (t, h)e^ (tf,to)x a . 

In particular, at t = tj, x^x(tf) = ||x a || 2 ^ 0, another contradiction. 

Thus in cither case we arrive at a contradiction, and so controllability implies 
the rank condition. 

Conversely, suppose that the system is not controllable. Then there exists an 
initial state xo £ R™ xl such that for all input signals u(t) £ E mxl , we have x(tf) ^ 
x f . Again, it follows from the commutativity of the convolution that 



Xf^x(tf) = / e A (t f ,a(s))Bu Xo (s) As + e A (tf,t )x 



to 



e A (s, t )Bu Xo (t f ,cr(s)) As + e A (tf,t )x 

3 

^ lk{s, t )A k Bu XQ (t f ,a(s)) As + e A (t f ,t )x - 
•» k=0 
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In particular, 



fc=0 



^^A k B I -f k (s,t )u Xo (t f ,a(s)) As ^ x f - e A (tf,t )x . 



Notice that the last equation holds if and only if there is no linear combination of 
the matrices A k B for k = 0, 1, . . . , n — 1, which satisfies 



n-l 

E 



A k Ba k = x f - e A (tf,t )x . 



The fact that there is no such linear combination follows from Lemma [2761 once we 
realize that an argument similar to the one given in the proof of this result holds if 
m < n. Thus, the matrix 

[B AB ■■■ A n ~ 1 B] 

cannot have rank n, and so we have shown that if the matrix has rank n, then it is 
controllable by contraposition. □ 

The preceding theorem is commonly called the Kalman Rank Condition after 
R. E. Kalman who first proved it in 1960 for the cases T = R and T = Z (see 
PS [20]). Therefore our analysis has unified the two cases, but we have also extended 
these results to an arbitrary time scale with bounded graininess. However, it is 
important to point out that the proof here is not the one that Kalman gave, which 
is the one classically used for K and Z (see [24] for example). In these two special 
cases, an observation about the particular form of the matrix exponential on R 
and Z (due to the uniform graininess) allows one to arrive at the result in a more 
straightforward manner. The general time scale case requires another argument 
altogether as demonstrated above. 

We now look at an example illustrating Theorem 12. 71 



Example 2.8. Consider the system 



- 8 


1 " 




2 






x{t) + 


1 


- 45 


10- 





u{t), x(0) = 
41 x(t). 



y(t) = [3 

It is straightforward to verify that 

rank [B AB] = rank ^ _i3/go 
so that the state equation is controllable by Theorem 12.71 



= 2, 



The next theorem establishes that there is a state variable change in the time 
invariant case that demonstrates the "controllable part" of the state equation. 

Theorem 2.9. Suppose the controllability matrix for the time invariant regressive 
linear system 

x A (t) = Ax{t) + Bu(t), x(t )=x , 
y(t)=Cx(t), 

satisfies 

rank [B AB ■■■ A^B] = q, 
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where < q < n. Then there exists an invertible matrix P such that 



P- Y AP = 



An 

(n-q)xq 



An 
M2 



P- X B = 



B n 

(n — q) xr 



where in is q x q, Bn is q x m, and 
rank [B u i u B u 



Al^Bn] = q- 



Proof. We begin constructing P by choosing q linearly independent columns pi, 
P2t ■ ■ , P q , from the controllability matrix for the system. Then choose p q +i, ■ ■ ■ ,p n 
as n x 1 vectors so that 

P=[Pl Pq Pq+1 Pn] 

is invertible. Define G so that PG = B. Writing the j-th column of B as a linear 
combination of the linearly independent columns of P given by pi,P2, ■ ■ ■ ,Pq, we 
find that the last n — q entries of the j-th column of G must be zero. This argument 
holds for j = 1, . . . , m, and so G = P~ 1 B does indeed have the desired form. 
Now set F = P~ 1 AP, yielding 

PF = [A Pl Ap 2 ■■■ A Pn ] . 

The column vectors Ap\ , . . . , Ap q can be written as linear combinations of p\ , . . . , p n 
since each column of A k B, k > can be written as a linear combination of these 
vectors. As for G above, the first q columns of F must have zeros as the last n — q 
entries. Thus, P _1 AP has the desired form. Multiply the rank-q controllability 
matrix by P^ 1 to obtain 

P- 1 [B AB ■■■ A n - 1 B] = [P- 1 B P- X AB ■■■ P^A^B] 

= [G FG ■■■ F^G] 

Bn AnBii ■ ■ ■ A\ Y 1 Bi 















Since the rank is preserved at each step, applying the Cayley-Hamilton theorem 
gives 

rank[B n A n B n ••• A q n 1 B 11 ]=q. 

a 

Next, we use the preceding theorem to prove the following. 
Theorem 2.10. The time invariant regressive linear system 
x A (t) = Ax(t) + Bu(t), x(t )=x , 
y(t)=Cx(t), 

is controllable if and only if for every scalar X the only complex n x 1 vector p 
satisfying p T A = Xp T , p T B = is p = 0. 

Proof. For necessity, note that if there exists p ^ and a complex A such that the 
equation given is satisfied, then 

p T [B AB ■■■ A n ~ 1 B] = [p T B p T AB ■ ■ ■ p T A n ~ 1 B~\ 

= [p T B Xp T B ■ ■ ■ X n - 1 p T B] , 

so that the n rows rows of the controllability matrix are linearly dependent, and 
hence the system is not controllable. 
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For sufficiency, suppose that the state equation is not controllable. Then by 
Theorem EH there exists an invertible P such that 



P-'AP : 



An 

0(n — <j) Xg 



Al2 
A 2 2 



P-'B 



P(n— q) xr 



with < q < n. Let p T = [Oi X q P^] P 1 , where p q is a left eigenvector for ^22 • 
Thus, for some complex scalar A, p^^4.22 = Xp^ , p q ^ 0. Then p ^ 0, and 



p T B = [0 p T q ] 







p T A = [0 p, 
Thus, the claim follows. 



An 




Ai 2 

A22. 



P-^[Q \p T q }P~ 



Xp 1 



□ 



The interpretation of Theorem 12.101 is that in a controllable time invariant sys- 
tem, A can have no left eigenvectors that are orthogonal to the columns of B. This 
fact can then be used to prove the next theorem. 

Theorem 2.11. The time invariant regressive linear system 

x A {t) = Ax(t) + Bu(t), x(t )=x , 

y{t) = Cx(t), 

is controllable if and only i/rank \zl — A B~\ = n for every complex scalar z. 

Proof. By Theorem 12.101 the state equation is not controllable if and only if there 
exists a nonzero complex n x 1 vector p and complex scalar A such that 

p T [XI- A B]p^ 0. 

But this condition is equivalent to rank [XI — A _B] < n. □ 

3. Observability 

Next, we turn our attention to observability of linear systems. As before, we 
treat the time varying case first followed by the time invariant case. 

3.1. Time Varying Case. In linear systems theory, when the term observability 
is used, it refers to the effect that the state vector has on the output of the state 
equation. As such, the concept is unchanged by considering simply the response of 
the system to zero input. Motivated by this, we define the following. 

Definition 3.1. The regressive linear system 

x A (t) = A(t)x(t), x(t ) = x , 
y(t) = C(t)x(t), 

is observable on [to,tf] if any initial state x(io) = %o is uniquely determined by the 
corresponding response y(t) for t G [to,tf). 

The notions of controllability and observability can be thought of as duals of 
one another, and so any theorem that we obtain for controllability should have an 
analogue in terms of observability. Thus, we begin by formulating observability in 
terms of an associated Gramian matrix. 
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Theorem 3.2 (Observability Gramian Condition). The regressive linear system 

x A (t) = A{t)x{t), x(t ) = x , 
y(t) = C(t)x(t), 

is observable on [to,tf] if and only if the n x n observability Gramian matrix 
Sfo(*o,t/) := f ' <i> T A (t,t )C T (t)C(t)<i> A (t,t )At, 

J to 

is invertible. 

Proof. If we multiply the solution expression 

y(t) = C(t)$ A {t,t )x , 
on both sides by ^^(t,t )C(t) and integrate, we obtain 

f ' $l(t,t )C T (t)y(t)At = &o(to,tf)x . 

J to 

The left side of this equation is determined by y(t) for t e [t ,tf), and thus this 
equation is a linear algebraic equation in xq. If ^o(to,tf) is invertible, then xq is 
uniquely determined. 

Conversely, if &o(to,tf) is not invertible, then there exists a nonzero vector x a 
such that &o{to,tf)x a — 0. But then x^^o(to,tf)x a = 0, so that 

C(t)$ A {t,t )x a = 0, t€[to,tf). 

Thus, x(t ) = x + x a yields the same zero-input response for the system as x(t ) = 
xq, and so the system is not observable on [to,i/]- □ 

The observability Gramian, like the controllability Gramian, is symmetric posi- 
tive semidcfinite. It is positive definite if and only if the state equation is observable. 

Once again we see that the Gramian condition is not very practical as it requires 
explicit knowledge of the transition matrix. Thus, we present a sufficient condition 
that is easier to check for observability. As before, observability and controllability 
can be considered dual notions to one another, and as such, proofs of corresponding 
results are often similar if not the same. Any missing observability proofs missing 
below simply indicates that the controllability analogue should be consulted. 

Definition 3.3. If T is a time scale such that n is sufficiently differentiable with the 
indicated rd-continuous derivatives, define the sequence of p x n matrix functions 

L (t) := C(t), 

L j+1 (t) := L^t)A{t) + Lf(t)(I + n(t)A(t)), j = 0, 1, 2, . . . 
As in the case of controllability, an induction argument shows that 

Lj(t) = ^-[cm A (t, S )]\ s=t . 

With this, an argument similar to the one before shows the following: 

Theorem 3.4 (Observability Rank Condition). Suppose q e Z + is such that, for 
t € [to,tf], C{t) is q-times rd- continuously differentiable and both of jj,(t) and A(t) 
are (q — 1) -times rd- continuously differentiable. Then the regressive linear system 

x A (t) = A{t)x{t), x(to)=x , 
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y(t) = C(t)x(t), 
is observable on [to,tf] if for some t c € [to, if), we have 



rank 



'Lo(t c j 
Li{t c ) 

L q {t c ) 



= n, 



wht 



L J (t) = ^-[cm A (t, S )] 



J = o, i, . 



3.2. Time Invariant Case. Like controllability, observability has equivalent con- 
ditions that become necessary and sufficient in the time invariant case. We begin 
with a Kalman rank condition for observability. 

Theorem 3.5 (Kalman Observability Rank Condition). The time invariant re- 
gressive linear system 

x A (t) = Ax(t) 1 x(t ) = x , 
y(t)=Cx(t), 

is observable on [to^f] if and only if the nm x n observability matrix 

C 

CA 



CA r 



has rank n. 



Proof. Again, we show that the rank condition fails if and only if the observability 
Gramian is not invertible. Thus, suppose that the rank condition fails. Then there 
exists a nonzero nxl vector x a such that 

CA k x a =0, k = 0,...,n-l. 

This implies, using Theorem 12.51 that 



&o(to,t f )x a = / el{t,t )C T Ce A (t,t )x a At 

Jtg 

e T A {t, t Q )C T Ik (t, t )CA k x a At 

3 k=0 

= 0, 

so that the Gramian is not invertible. 

Conversely, suppose that the Gramian is not invertible. Then there exists 
nonzero x a such that xj^o(to, tf) x a — 0. As argued before, this implies 

Ce A (t, t a )x a = 0, t€[t ,tf). 
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At t = to, we obtain Cx a — 0, and differentiating k times and evaluating the resuit 
at t = to gives CA k x a = 0, k = 0, . . . , n — 1. Thus, 

C 



CA 



0, 



and the rank condition fails. 



□ 



The proof of the preceding result demonstrates an important point about con- 
trollability and observability in the arbitrary time scale setting: namely, proofs 
of similar results for the two notions are often similar, but can sometimes be very 
different. Indeed, comparing the proof of the Kalman condition for controllability 
with the proof of the Kalman condition for observability highlights this contrast. 

The following example uses Theorem 13.51 



Example 3.6. Consider the system 
x A (t) = 



- 8 


1 " 




2" 


45 


3 °1 


x(t) + 


1 


. 45 


10. 





*(*), z(0) = 



y(t)= [3 i]x(t). 

From Example 12.81 recall that the system is controllable. We claim the system is 
also observable. This follows from Theorem 13.51 since 

C 



rank 



CA 



= rank 



28 
45 



The following three theorems concerning observability have proofs that mirror 
their controllability counterparts, and so will not be given here. 

Theorem 3.7. Suppose the observability matrix for the time invariant regressive 
linear system 

x A (t) = Ax{t) + Bu(t), x(t )=x , 
y(t)=Cx(t), 



satisfies 



rank 



C 
CA 



CQ = [C u 0] 



CA 71 - 1 

where < £ < n. Then there exists an invertible n x n matrix Q such that 

A.21 A22_ 

where An is £ x i, C\\ is p x I, and 

Cu 
CiiAn 



rank 



= £. 
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The state variable change Theorem 13.71 is constructed by choosing n — I vectors 
in the nullspace of the observability matrix, and preceding them by i vectors that 
yield a set of n linearly independent vectors. 

Theorem 3.8. The time invariant regressive linear system 
x A (t) = Ax(t) + Bu(t), x(t )=x , 
y(t) = Cx(t), 

is observable if and only if for every complex scalar X, the only complex n x 1 vector 
p that satisfies Ap = Xp, Cp = is p = 0. 

Again, Theorem l3.8l can be restated as saying that in an observable time invariant 
system, A can have no right eigenvectors that are orthogonal to the rows of C. 

Theorem 3.9. The time invariant regressive linear system 
x A {t) = Ax(t) + Bu{t), x(t )=xo, 
y(t) = Cx(t), 

is observable if and only if 

= n, 

for every complex scalar z. 

4. Realizability 

In linear systems theory, the term realizability refers to the ability to characterize 
a known output in terms of a linear system with some input. We now make this 
precise. 

Definition 4.1. The regressive linear system 

x A = A(t)x{t) +B(t)u(t), x(t o ) = 0, 
y(t) = C(t)x(t), 

of dimension n is a realization of the weighting pattern G(t, cr(s)) if 

G(t,a{s)) =C(t)$ A (t,o-(s))B{s), 

for all t, s. If a realization of this system exists, then the weighting pattern is 
realizable. The system is a minimal realization if no realization of G(t,a(s)) with 
dimension less than n exists. 

Notice that for the system 

x A (t) = A{t)x(t) + B{t)u{t), x(t Q ) = 0, 

y(t)=C(t)x(t)+D(t)u(t), 

the output signal y(t) corresponding to a given input u(t) and weighting pattern 
G{t,o-(s)) = C(t)$ A {t,a(s))B{s) is given by 

y(t)= / G(t,a(s))u(s)As + D(t)u(t), t>t . 

J t 

When there exists a realization of a particular weighting response G(t, cr(s), there 
will in fact exist many since a change of state variables will leave the weighting 
pattern unchanged. Also, there can be many different realizations of the same 



rank 



C 

zl - A 
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weighting pattern that all have different dimensions. This is why we are careful to 
distinguish between realizations and minimal realizations in our definition. 

We now give equivalent conditions for realizability: as before, we begin with the 
time variant case and then proceed to the time invariant case. 

4.1. Time Varying Case. The next theorem gives a characterization of realizable 
systems in general. 

Theorem 4.2 (Factorization of G(t,a(s))). The weighting pattern G(t,a(s)) is 
realizable if and only if there exist a rd-continuous matrix H(t) that is of dimension 
q x n and a rd-continuous matrix F(t) of dimension n x r such that 

G(t, a(s)) = H(t)F(a(s)), for all t, s. 

Proof. Suppose there exist matrices H(t) and F(t) with G(t 7 cr(s)) = H(t)F(a(s)). 
Then the system 

x A (t) =F(t)u(t), 
y(t) = H(t)x{t), 

is a realization of G(t, a(s)) since the transition matrix of the zero system is the 
identity. 

Conversely, suppose that G(t, a(s)) is realizable. We may assume that the system 
x A {t) = A(t)x(t) + B(t)u(t), 
y(t) = C(t)x(t), 

is one such realization. Since the system is regressive, we may write 

G{t,a(s)) = C{t)$ A (t,o-(s))B(s) = C(t)$ A (t,0)$A(0,o-(s))B(s). 
Choosing H(t) = C(t)$ A (t, 0) and F(t) = $ A (0, a(t))B(t), the claim follows. □ 

Although the preceding theorem gives a basic condition for realization of linear 
systems, often in practice it is not very useful because writing the weighting pattern 
in its factored form can be very difficult. Also, as the next example demonstrates, 
the realization given by the factored form can often be undesirable for certain 
analyses. 

Example 4.3. Suppose T is a time scale with < \i < 2. Under this assumption, 
-1/4 e U + (T). Then the weighting pattern 

G(t,a(s)) = e_i /4 (t,<7(s)), 

has the factorization 

G(t,a(s)) = e_ 1/4 (i,cr(s)) = e_ 1/4 (t, 0)e e( _ 1/4) (cr(s), 0). 

By the previous theorem, a one-dimensional realization of G is 

x A (t) = e e( _ 1/4) (t,0)w(i), 

y(t) = e_ 1/4 (t,0)x{t). 

This state equation has an unbounded coefficient and is not uniformly exponentially 
stable (note that e e (_i/ 4 )(f, 0) = ei/( 4 _^)(t, 0) is unbounded since 1/(4 — y) > 0). 
However, the one-dimensional realization of G given by 

x A (t) = ~x(t)+u(t), 
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y(t) = x(t), 

does have bounded coefficients and is uniformly exponentially stable. 

Before examining minimal realizations, some remarks are in order. First, note 
that the inverse and a operators commute: 

= p-\t) + ^t)(-P(a(t))r 1 P A (t)P-\t) 
= P-\t) (P(a(t))-\P(a(t)) Pi^P- 1 ® 
= (P(a(t)))-K 

Second, it is possible to do a variable change on the system 
x A (t) =A(t)x(t) +B(t)x(t), 
y(t) = C(t)x(t), 

so that the coefficient of x A (t) in the new system is zero, while at the same time 
preserving realizability of the system under the change of variables. 
Indeed, set z{t) = P~ 1 (t)x(t) and note that P(t) = &A(t,to) satisfies 

(P(a(t)))-'A(t)P(t) (P(a(t)))- 1 P A (t) = 0. 

If we make this substitution, then the system becomes 

z A (t) = p-\a(t))B(t)u(t), 

y(t) = C(t)P{t)z(t). 

Thus, in terms of realizability, we may assume without loss of generality that A{t) = 
by changing the system to the form given above. 

It is important to know when a given realization is minimal. The following 
theorem gives a necessary and sufficient condition for this in terms of controllability 
and observability. 

Theorem 4.4 (Characterization of Minimal Realizations). Suppose the regressive 
linear system 

x A (t) =A(t)x(t) + B(t)x(t), 
y{t) = C(t)x(t), 

is a realization of the weighting pattern G(t, cr(s)). Then this realization is minimal 
if and only if for some t and tf > t the state equation is both controllable and 
observable on [t ,tf]. 

Proof. As argued above, we may assume without loss of generality that A(t) = 0. 
Suppose the n-dimensional realization given is not minimal. Then there is a lower 
dimension realization of G(t, cr(s)) of the form 

z A (t) = R(t)u(t), 
y(t) - S{t)z(t), 

where z{t) has dimension n z < n. Writing the weighting pattern in terms of both 
realizations produces C(t)B(s) = S(t)R(s) for all t, s. Thus, 

C T (t)C(t)B(s)B T (s) = C T (t)S(t)R(s)B T (s), 
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for all t, s. For any to and tf > to, it is possible to integrate this expression with 
respect to t and then with respect to s to obtain 



&o(to,tf)&c(to,t f ) = f ' C T (t)S(t)At [ ' R(s)B T {s)As. 

Jt Jt a 

The right hand side of this equation is the product of an n x n z matrix and an 
n z xn matrix, and as such, it cannot have full rank since the dimension of the space 
spanned by the product is at most n z < n. Therefore, ^o(to,tf) and ^c{to,tf) 
cannot be simultaneously invertible. The argument is independent of the to and tf 
chosen, and so sufficiency is established. 

Conversely, suppose that the given state equation is a minimal realization of the 
weighting pattern G(t, c(s)), with A(t) = 0. We begin by showing that if either 

Sfc(*o,t/) = / f B(t)B T (t) At, 



to 



or 



tf 



&o(to,t f )= / C 1 (i)C(t) At, 

Jt 

is singular for all to and tf, then minimality is violated. Thus, there exist intervals 
[to,t a f ] and [to,t b f ] such that <^ c (to,t a f ) and & (to,t b f ) are both invertible. If we 
let t = mm {ioS*o} and tf — max{t^,iy}, then the positive definiteness of the 
observability and controllability Gramians yields that both ^c{toi tf) and ^o{to,tf) 
are invertible. 

To show this, we begin by supposing that for every interval [to,tf] the matrix 
&c(to,tf) is not invertible. Then, given to and tf there exists an n x 1 vector 
x = x(to, tf) such that 

ftf 

= x Tc S c (to,tf)x = / B{t)B T {t)x At. 

J t 

Thus, x T B(t) = for t £ [t ,t/). 

We claim that there exists at least one such vector x that is independent of to 
and tf. To this end, note that if T is unbounded from above and below, then for 
each positive integer k there exists annxl vector Xk with 

||x fc || = 1, xlB(t)=0, t£[-k,k]. 

Thus, {xk}^Li is a bounded sequence of n x 1 vectors and by the Bolzano- Wicr- 
strauss Theorem, it has a convergent subsequence since T is closed. We label this 
convergent subsequence by {xkj }j X Li and denote its limit by x = lim^oo x/^ . Note 
x^B(t) — for all t, since for any given time t a , there exists a positive integer J a 
such that t a £ [—kj, kj] for all j > J a , which in turn implies x^._B(t a ) = for all 

j > J a - Hence, Xq satisfies x^B(t a ) = 0. 

Now let be a constant, invertible, n x n matrix with bottom row Xq. Using 
P^ 1 as a change of state variables gives another minimal realization of the weighting 
pattern, with coefficient matrices 



P^Bit) = 



Bi(t) 

Olxm 



C(t)P=[C 1 (t) C 2 (t)], 
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where Bi(t) is (n — 1) X m, and C\(t) is p x (n — 1). Then a straightforward 
calculation shows G(t,a(s)) = C\{t)Bi{a{s)) so that the linear state equation 

z A (t)=B 1 {t)u(t), 
y{t) = C 1 {t)z{t), 

is a realization for G(t, er(s)) of dimension n — 1. This contradicts the minimality 
of the original n-dimensional realization. Thus, there must be at least one ip and 
one tf > *q such that ^c?(*o>*/) is invertible. 

A similar argument shows that there exists at least one t\ and one tf > t\ 
such that &o(t\,t b f) is invertible. Taking to = rninjtQ, fg} and f/ = max{ty,ty} 
shows that the minimal realization of the state equation is both controllable and 
observable on [to,*/]. □ 

4.2. Time Invariant Case. We now restrict ourselves to the time invariant case 
and use a Laplace transform approach to establish our results. Instead of consid- 
ering the time-domain description of the input-output behavior given by 

y(t)= f G(t,a(s))u(s)As, 
Jo 

we examine the corresponding behavior in the z-domain. Laplace transforming the 
equation above and using the Convolution Theorem [T3] yields Y(z) = G(z)U(z). 
The question is: given a transfer function G(z), when does there exist a time 
invariant form of the state equation such that 

C(zl- A)- 1 B = G(z), 

and when is this realization minimal? 

To answer this, we begin by characterizing time invariant realizations. In what 
follows, a strictly-proper rational function of z is a rational function of z such that 
the degree of the numerator is strictly less than the degree of the denominator. 

Theorem 4.5. Thepx q transfer function G(z) admits a time invariant realization 
of the regressive linear system 

x A {t) = Ax(t)+Bu(t), 
y(t)=Cx(t), 

if and only if each entry of G{z) is a strictly-proper rational function of z. 

Proof. If G(z) has a time invariant realization, then G has the form G(z) = C(zl — 
A)~ 1 B. We showed in [TJJ [TJ] that for each Laplace transformable function /(t), 
F(z) — > as z — > co, which in turn implies that F(z) is a strictly-proper rational 
function in z. Thus, the matrix (zl — A)~ l is a matrix of strictly-proper rational 
functions, and G{z) is a matrix of strictly-proper rational functions since linear 
combinations of such functions are themselves strictly-proper and rational. 

Conversely, suppose that each entry Gij(z) in the matrix is strictly-proper and 
rational. Without loss of generality, we can assume that each polynomial in the 
denominator is monic (i.e. has leading coefficient of 1). Suppose 

d(z) = z r + d r ^ 1 z r ^ 1 H h d 
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is the least common multiple of the polynomials in denominators. Then d(z)G(z) 
can be decomposed as a polynomial in z with p x q constant coefficient matrices, 
so that 



d(z)G(z) = P r -iz r - x + ■ ■ ■ + P X Z + Pq 



We claim that the qr- 


dimensional matrices 


given by 






' % 


h 


% 








% 


% 


% 




% 


A = 








, B = 






% 


% 


% 




% 




-doI q 


-d\I q 






A 


form a realization of G(z). 


To see this, let 












R{z) = (zl- 


A)~ l B, 





c 



Po 
Pi 

P, — i 



and partition the qr x q matrix R(z) into r blocks Ri(z), R%{z), . . . ,R r (z), each 
of size q x q. Multiplying R(z) by [zl — A) and writing the result in terms of 
submatrices gives rise to the relations 

Ri+i(z) = zRi, i = l,...,r — 1, 

zR r (z) + d Ri(z) + diR 2 {z) H h d r - 1 R r (z) = I q . 

Using (14. ip to rewrite (14.21) in terms of R\(z) gives 



(4.1) 
(4.2) 



Ri{z) 

and thus from (|4.ip again, we have 



R(z) 



1 



d{z) Iq ' 



Zl„ 



d(z) 



Multiplying by C yields 

CizI - Ay x B 
which is a realization of G(z). 



d(z) 



{Po + zPi 



l P r . x ) = G(z), 



□ 



The realizations that are minimal are characterized in the following theorem, 
which is repeated here for completeness sake in this easier case of Theorem 14.41 

Theorem 4.6. Suppose the time invariant regressive linear system 

x A (t) = Ax(t)+Bx(t), 
y(t)=Cx(t), 

is a realization of the transfer function G(z). Then this state equation is a minimal 
realization of G(z) if and only if it is both controllable and observable. 

Proof. Suppose the state equation is a realization of G(z) that is not minimal. Then 
there is a realization of G(z) given by 

z A {t) = Pz(t) + Qz(t), 

y(t) = Rz(t), 
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with dimension n z < n. Thus, 

Ce A (t, 0)B = Re P (t, 0)Q, t > 0. 
Repeated differentiation with respect to t, followed by evaluation at t = yields 

CA k B = RF k Q, fc = 0,l,... 
Rewriting this information in matrix form for k = 0, 1, . . . , In — 2, we see 
CB CAB ■■■ CA n ~ l B 



_CA n - x B CA n B ••• CA 2n - 2 B 
which can be rewritten as 
C 

1 

A n ~ 1 B 



RQ RPQ 
RP n ~ 1 Q RP n Q 



RP n - x Q' 

R p2n-2Q 



CA 



CA 



n-l 



[B AB ■■ 



R 

RP 



RP 



■n— 1 



[Q PQ 



yn — 1 



Q] ■ 



However, since the right hand side of the equation is the product of an n z p x n z 
and an n z x n z m matrix, the rank of the product can be no greater than n z . 
Thus, n z < n, which implies that that the realization given in the statement of the 
theorem cannot be both controllable and observable. Therefore, by contraposition 
a controllable and observable realization must be minimal. 

Conversely, suppose the state equation given in the statement of the theorem is 
a minimal realization that is not controllable. Then there exists an n x 1 vector 
y 7^ such that 

y T [B AB ■■■ A^B] = 0, 

which implies y T A k B = for all k > by the Cay ley-Hamilton theorem. For 
P^ 1 an invertible n x n matrix with bottom row y T , then a variable change of 
z{t) = P~ 1 x(t) produces the state equations 

z A (t) = Az(t) + Bu(t), 
y(t) - Cz(t), 

which is also an n-dimensional minimal realization of G(z). Partition the coefficient 
matrices of the state equation above as 



A = P^AP 



An A 12 
A21 A 2 2 



B = R- Y B 



B x 




C = CP = [Ct CA] , 



where An is (n — 1) x (n — 1), B\ is (n — 1) x 1, and C\ is 1 x (n — 1). From these 
partitions, it follows from the construction of P that AB — P~ 1 AB has the form 



'AnBj 




'AnB{ 


A 2 \Bi_ 








AB = 

Since the bottom row of P~ 1 A k B is zero for all k > 



A k B = 



-a\ x b; 





k > 0. 
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But, An, Bi,C± give an (n — l)-dimcnsional realization of G(z) since 



Ce A (t,0)B = [Ci C 2 ]Y j A k Bh k {t,Q) 



k=0 



Co 



k=0 



= Cie A Jt,Q)B x 



A\ x Bi 




h k {t, 0) 



so that the state equation in the statement of the theorem is in fact not minimal, 
a contradiction. A similar argument holds if the system is assumed not to be 
observable. □ 



We now illustrate Theorem 14.51 and Theorem 14.61 with an example. 
Example 4.7. Consider the transfer function 

_ 9(37 + 300z) 
[Z) ~ 5 + 75z + 270z 2 ' 
G{z) admits a time invariant realization by Theorem \4-5\ since G{z) is a strictly- 
proper rational function of z. The form of G(z) indicates that we should look for a 
2-dimensional realization with a single input and single output. We can write 

-l 



G{z) = [3 41 zl 



1 

45 



10 



so that a time invariant realization of G(z) is given by 



x A (t) 



- 8 


1 " 




2 




30 1 


x{t) + 


1 


. 45 


10_ 





u(t), x(0) = Xq, 



y(t)=[3 4]z(i). 



We showed in Examvle \2. 8\ that this realization is in fact controllable, and we showed 
in Example \3.6\ that it is also observable. Thus, Theorem \4-6\ guarantees that this 
realization of G(z) is minimal. 

5. Stability 

We complete our foray into linear systems theory by considering stability. In 
|23j . Potzsche, Siegmund, and Wirth deal with exponential stability. DaCunha 
also deals with this concept under a different definition in [TOl [TT] and emphasizes 
the time varying case. We begin by revisiting exponential stability in the time 
invariant case and then proceed to another notion of stability commonly used in 
linear systems theory. 

5.1. Exponential Stability in the Time Invariant Case. We start this section 
by revisiting the notion of exponential stability. We are interested in both the 
time invariant and time varying cases separately since it is often possible to obtain 
stronger results in the time invariant case. 

We have already noted that if A is constant, then $^(t,to) = e^^, to). In what 
follows, we will consider time invariant systems with <o = in order to talk about 
the Laplace transform. 

DaCunha defines uniform exponential stability as follows. 
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Definition 5.1. [II] The regressive linear system 

x A = A(t)x(t), x(t ) = Xq, 

is uniformly exponentially stable if there exist constants 7, A > with —A G 1Z + 
such that for any to and x(io), the corresponding solution satisfies 

Mt)\\ < ||aj(t )||7e-A(t,*o), t>t . 

With this definition of exponential stability, we can prove the next theorem. 

Theorem 5.2. The time invariant regressive linear system 

x A (t) = Ax(t), x(to) = xq, 

is uniformly exponentially stable if and only if for some /3 > 0, 

\\e A (t,t a )\\At < (3. 

to 

Proof. For necessity, note that if the system is uniformly exponentially stable, then 
by [H3 Theorem 3.2], 



f ||e A (t,0)||At< f 
Ja Jo 



7 e_ A (t,0)At = J, 



so that the claim follows. 

For sufficiency, assume the integral condition holds but for the sake of contra- 
diction that the system is not exponentially stable. Then, again by [TOl Theorem 
3.2], for all A, 7 > with -A S TZ+ , we have ||e A (t,0)|| > 7e_ A (i,0). Hence, 



' (Mil A/ • / 7 e_ A (t,0)At = X e _ A ( t) o) 



2 
A' 



In particular, if we choose 7 > /3A, then 

^°°||e A (t ) 0)||At>^=^ 
a contradiction. □ 
Now consider the system 

x A (t) = Ax(t), x(0) = I. 
Transforming this system yields 

X{z) = {zI-A)-\ 

which is the transform of eA{t, 0). This result is unique as argued in [T31IIH]- Note 
that this matrix contains only strictly-proper rational functions of z since we have 
the formula 

v ' det(z/ - A) 

Specifically, det(z/ — A) is an nth degree polynomial in z, while each entry of 
adj(z/ — ^4) is a polynomial of degree at most n — 1. Suppose 

det(zl -A) = (z- AO^ 1 (z - \ 2 )^ • • • (z - A m )^, 
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where Ai, A2, . . . , X n are the distinct eigenvalues of the n x n matrix A, with corre- 
sponding multiplicities ipi, ip2, ■ ■ ■ , ipm- Decomposing (zl — A) -1 in terms of partial 
fractions gives 

m V>fc 

fc=ij=i v ; 

where each Wkj is an n x n matrix of partial fraction expansion coefficients given 

by 

If we now take the inverse Laplace transform of (zl — A)^ 1 in the form given above, 
we obtain the representation 

o) = E E %%rr e ^ (*• °). ( 5 - x ) 

where fj(fi, \k) is the sequence of functions obtained from the residue calculations 
of the jth derivative in the inversion formula. For example, the first few terms in 
the sequence are 



fo(pAk) = 1, 



4 1 



7 1 + fiX k 



(f — 

Wo 1 + P 



2 



/ 2 (M*)=( / ^^-Ar) - / ^ / x At, 



, — — v) Ar /* 

Mfc ' Jo (l + ^Afe) Jo 1 + Mfc 



At, 



(1 + Mfc) 



Notice that if /1 is bounded, then each fj(fi, Afe) can be bounded by a "regular" 
polynomial of degree j in t, call it a,j(t). That is, /j can be bounded by functions 
of the form cij(t) = cijP + aj-iP ■ ■ ■ + clq. This observation will play a key role 
in the next theorem. Potzsche, Siegmund, and Wirth do prove this result in [23], 
but our proof differs from theirs in that we use new transform results to obtain 
it, while they use other techniques. Note, however, that in the theorem we do use 
their definition of exponential stability rather than the one given by DaCunha. For 
completeness, we remind the reader by restating their definition here. 

Definition 5.3. [23] For t,t e T and a; G K™, the system 

x A = A(t)x, x(to) = xq, 

is uniformly exponentially stable if there exists a constant a > such that for every 
to ET there exists a K > 1 with 

HMMo)|| < JCe-^*-*") fori>i , 
with K being chosen independently of to. 
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Recall that DaCunha's definition of uniform exponential stability of a system will 
imply that the system is uniformly exponential stable if we use Potzsche, Siegmund, 
and Wirth's definition of the concept, but the converse need not be true in general. 
Thus, DaCunha's definition is weaker in this sense. 

Theorem 5.4 (Spectral Characterization of Exponential Stability). Let T be a 

time scale which is unbounded above but has bounded graininess. The time invariant 
regressive linear system 

x A (t) = Ax(t) 1 x(t ) = x , 

is uniformly exponentially stable (in the sense of Definition \ 5.3$ if and only if 
spec(^4) C <S(T), the regressive set of exponential stability for T, given by 



S(T) := {A e 



1 



lim sup 

T^oo 2 — to J to s\fi(t) 



' lim lQg ' 1 + sA 'At<0 } . 



Proof. Suppose the eigenvalue condition holds. Then, appealing to Theorem 
and writing the exponential in the explicit form given above in terms of the distinct 
eigenvalues Ai, A2, • ■ • , A m , we obtain 



|e A (t,0)||At = 







rn ipk r I % 

EE^%%#^ <'- ( » 

^=1 0=1 



m ipk 

< EE 

k=i 0=1 

m ipk 

< EE 11^ 

fe=i j=i 

m ipk 



(3-1)1 
fj-i((i, Xk] 



At 



< EE / 

m ipk „ c 

< EE 11^11 / 



(i-i)! 

|o i _ 1 (t)e Alk (t,0)| At 
aj _i(t)e- Qt At 



At 



ij-_i(t)e a * dt < 00. 

Thus, by Theo- 



fc=i 3=1 

Note that the last three lines hold by appealing to Definition! 
rem 15.21 the system is uniformly exponentially stable. 

Now, for the sake of a contradiction, assume that the eigenvalue condition fails. 
Let A be an eigenvalue of A with associated eigenvector v, with A ^ 5(C). Direct 
calculation shows that the solution of the system x A — Ax, x(0) — v, is given by 
x(t) = e\(t, Q)v. From [23] . if A ^ 5(C), then lim^oo e\(t, 0) ^ 0, so that we arrive 
at a contradiction. □ 

5.2. BIBO Stability in the Time Varying Case. Besides exponential stability, 
the concept of bounded-input, bounded-output stability is also a useful property 
for a system to have. As its name suggests, the notion is one that compares the 
supremum of the output signal with the supremum of the input signal. Thus, we 
define the term as follows. 

Definition 5.5. The regressive linear system 

x A (t) = A(t)x{t) + B(t)u(t), x{t ) = x , 
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y(t) = C(t)x(t), 

is said to be uniformly bounded-input, bounded- output (BIBO) stable if there exists 
a finite constant i] such that for any to and any input u(t) the corresponding zero- 
state response satisfies 

sup \\y(t)\\ < r/sup ||u(i)||. 

t>t t>t 

Note that we use the word "uniform" to stress that the same n works for all to 
and all input signals. 

The following characterization of BIBO stability is useful. 

Theorem 5.6. The regressive linear system 

x A (t) = A(t)x(t) + B{t)u(t), x(to) = x , 
y(t) = C(t)x(t), 

is uniformly bounded-input, bounded- output stable if and only if there exists a finite 
constant p such that for all t, r with t>r, 

' \\G(t,a(s))\\As<p. 



i 



[ C(t)$ A (t,o-(s))B(s)u(s)As <f ||G(t,<r(«))||||u(s)||As, t>t Q 

J to J tQ 



Proof. Assume such a p exists. Then for any to and any input signal, the corre- 
sponding zero-state response of the state equation satisfies 

ll</(/>|| = 

Replacing by its supremum over s > t , and using the integral condition, 

we obtain 

\\y(t)\\ < sup||u(t)|| / \\G(t,a(s))\\ As < psup ||«(t)||, t > t . 

t>t Jt t>t 

Thus, the system is BIBO stable. 

Conversely, suppose the state equation is uniformly BIBO stable. Then there 
exists a constant r\ so that, in particular, the zero-state response for any to and any 
input signal such that sup t>to ||u(f)|| < 1 satisfies sup t>to ||y(f)|| < i]. For the sake 
of a contradiction, suppose no finite p exists that satisfies the integral condition. 
Then for any given p > 0, there exist r p and t p > t p such that 

"\\G(t p ,a(s))\\As>p. 

p 

In particular, if p = n, this implies that there exist t v , with t v > t v , and indices 
i, j such that the i, j-entry of the impulse response satisfies 

\G l3 (t v ,a(s))\ As > T]. 

With to = t v consider the m x 1 input signal u(t) defined for t > to as follows: set 
u(t) = for t > t v , and for t G [^O;^] set every component of u(t) to zero except 
for the j-th component given by the piecewise continuous signal 

fl, G ij (t ri ,a(t))>0, 
Uj (t):=lo, Gij(t v , a(t)) = 0, 
1-1, G tf (t JJ ,tr(t))<0. 
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This input signal satisfies ||u(t)|| < 1 for all t > to, but because of the integral con- 
dition above, the i-th component of the corresponding zero-state response satisfies 

yi(t v )= G i j(t r „a(s))u j (s)As = \Gij(t v ,a(s))\As>r]. 

J to J to 

Since 112/(^)11 > we arrive at a contradiction that completes the proof. □ 

We now wish to give conditions under which the notions of exponential stability 
and BIBO stability are equivalent. To this end, we begin with the following. 

Theorem 5.7. Suppose the regressive linear system 

x A (t) = A(t)x{t) + B(t)u(t), x(t ) = 350, 
y(t) = C(t)x(t), 

is uniformly exponentially stable, and there exist constants (3 and 7 such that 
11-5(011 — an d l|C(*)ll < a f or a M t- Then the state equation is also uniformly 
bounded-input, bounded- output stable. 

Proof. Using the bound implied by uniform exponential stability, we have 
\G(t, a(s))\\ As < J* \\C(t)\\ \\<t> A (t, a(s))\\ \\B(s)\\ As 

<a(3 f \\<S> A (t,a(s))\\As 



<a/3 J 7e_ A (i,cr(s))As 
a/? 7 /"* A 



< 



A J T 1 — n(s)\ 

a/?7 



x/(i-fiX)(s,t) As 



N (l-e_ A (t,r)) 
< a/37 



A 

By Theorem 15. 61 the state equation is also bounded-input, bounded-output stable. 

□ 

The following example illustrates the use of Theorem 1 5. 71 
Example 5.8. Let T be a time scale with < \i < |. Consider the system 



x A (t) 



-2 1 
-1 -sin(i) 



x(t) 



cos(i) 
sin(t) 



u(t), x(t )=x , 



y(t)=[l e^(t,0)]x(t), 



where here, sin(t) and cos(i) are the usual trigonometric functions and not their 
time scale counterparts. DaCunha shows that the system is uniformly exponentially 
stable by applying [TU1 Theorem 4.2] (also found in [TTJ Theorem 3.2]) with the 
choice Q(t) = I. For t > 0, we have \\B(t)\\ = 1 and ||C(f)|| = V 1 + (e_i(i,0)) 2 < 
\J~2, since p = —1 G 1Z + from our assumption on T. By Theorem 15 .7\ the state 
equation is also uniformly bounded-input, bounded-output stable. 



28 J. M. DAVIS, I. A. GRAVAGNE, B. J. JACKSON, R. J. MARKS II 



EJDE-2009/37 



For the converse of the previous theorem, it is known on T = E and T = Z that 
stronger hypotheses than simply having the system be BIBO stable are necessary to 
establish exponential stability (see [TJ [^J, [221 EI])- At present, we lack an analogue 
of this result for an arbitrary time scale in the time varying system case. We will 
see that the time invariant case does allow for the equivalence of the two notions 
in the general time scale case under certain conditions. 

5.3. BIBO Stability in the Time Invariant Case. In order to extend the 
definition of BIBO stability to the time invariant case, we first need the following 
definitions. 

Definition 5.9. [13l[19] Let u(t) € C prd -o2(T, M). The shift ofu(t) by a(s), denoted 
by u(t,a(s)), is given by 

u(t,O=J?- 1 {U(z)e° ez (s,0)}, 

where U(z) :— ££ {u(t)}(z), and J?f, J?f -1 denote the generalized Laplace transform 
and its inverse. 

Definition 5.10. [T3l [T9] For /, g 6 C pr d-c2 (T, K) , the convolution f* g is given by 

(/*<?)(*)= / f(r)g(t,a(T))AT. 
Jo 

Definition 5.11. For any shift u(t,a(s)) of the transformable function u(t), the 
time invariant system 

x A (t) = Ax(t) +Bu(t), x(t )=x , 
y(t)=Cx(t), 

is said to be uniformly bounded-input, bounded- output stable if there exists a finite 
constant n such that the corresponding zero-state response satisfies 

sup||y(t)|| < ?7supsup || u(t,a(s)) || . 

t>0 t>0 s>0 

Note that Definitions 15 . 51 and 1 5 . 1 1 1 are different: one deals with the time varying 
case and the other with the time invariant case. The modified definition in the time 
invariant case says that the output stays bounded over all shifts of the input. 

Theorem 5.12. The time invariant regressive linear system 

x A (t) = Ax{t) + Bu(t), x(t )=x , 
y(t) = Cx(t), 

is bounded-input, bounded- output stable if and only if there exists a finite (3 > 
such that 

\\G{t)\\ At < f3. 

Proof. Suppose the claimed f3 > exists. For any time t, we have 

y(t)= [ Ce A (t,a(s))Bu(s) As = [ Ce A (s, 0)Bu(t, cr(s)) As, 
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h(t)\\ < \\C\\ ||e x ( a ,0)||||B|| snp \\u(t, cr(s))\\ As 

0<s<t 



29 



< \\c\> 



M«,0)||As||B|| sup ||u(i, a(s))\\, 

8>0 



which implies 

roc 

sup\\y(t)\\ < \\C\\ / ||eA(*,0)||A«||B|| supsup ||u(i, a(s))\\. 

t>0 Jo t>0 s>0 

If we choose r? = ||C|| (3 \\B\\, the claim follows. 

Conversely, suppose that the system is bounded-input bounded-output stable, 
but for the sake of a contradiction that the integral is unbounded. Then, 



sup||y(»|| < 77supsup||w(i,cr(s))||, 

t>0 t>0 s>0 



and 



f 

Jo 



\\G(t)\\At > 13 for all (3 > 0. 



In particular, there exist indices i,j such that 

|Gy(t)|At > 0. 



I 

Jo 



Choose u(t, <t(s)) in the following manner: set Uk(t,a(s)) = for all k ^ j, and 
define Uj(t,a(s)) by 

fl, ifGy(«)>0, 

Uj(t,a{s)) := lo, ifG y («)=0, 

[-1, ifGy(«)<0. 

Choose (3 > 7] > 0. Note sup t > sup s > \\u(t, a(s)\\ < 1, so sup t > \\y(t)\\ < ij. 
However, 



sxi-p\\y( 

t>o 



sup 



G(s)u(t,a(s))As 



= sup 

t>o 



Gj(s) ■ Uj(s) As 

Jo 

/ \Gij{8)\As 
Jo 



> sup 

t>0 



>/3>?7, 

which is a contradiction. Thus, the claim follows. 



□ 



The next theorem demonstrates the equivalence of exponential and BIBO sta- 
bility in the time invariant case. Recall that this is a notion we currently lack in 
the time varying case. 



.'ill 
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Theorem 5.13 (Equivalence of BIBO and Exponential Stability). Suppose the 
time invariant regressive linear system 

x A (t) = Ax(t) + Bu(t), x(t )=x , 

y(t) = Cx(t), 

is controllable and observable. Then the system is uniformly bounded-input, bounded 
output stable if and only if it is exponentially stable. 

Proof. If the system is exponentially stable, then 



poo poo 

/ \\Ce A (t,0)B\\At<\\C\\\\B\\ ||e A (t, 0)|| At 
Jo Jo 



< 



V, 



by Theorem 15.21 

Conversely, suppose the system is uniformly bounded-input, bounded output 
stable. Then 



which implies 



/ \\Ce A (t,Q)B\\ At < oo, 
Jo 



lim Ce A {t,0)B = 0. 

t— >GO 



(5.2) 



(5.3) 



Using the representation of the matrix exponential from ()5.1|) , we may write 

Ce A M)B = £ £ N *J (j-lV ^ (*' °) > 

where the Xk are the distinct eigenvalues of A, the Nkj are constant matrices, and 
the fj(n, Xk) are the terms from the residue calculations. In this form, 



— Ce A (t,0)B 



' fi(^A ( i!! /'I/') A/, I A,/, ,(/i.A, ) 

fc=l ' j=2 



= EK+E( : 



(7-2)! 



(7-1)! 



ex k (t,0). 



If this function does not tend to zero as t — > oo, then using (|5.3[) and (|5.2p , we 
arrive at a contradiction. Thus, 

lim (-^-Ce^M)^ = lim Cy4.e A («, 0)B = lim Ce A (t, 0)AB = 0, 

t — >oo V / t — >oo t — too 

where the last equality holds by noting that if A is constant, then A and e A (t,0) 
commute. Similarly, it can be shown that any order time derivative of the expo- 
nential tends to zero as t — > oo. Thus. 



lim CA t e A (t,0)A J B = 0, i,j = 0, 1, 



It follows that 



lim 

t — ► oc 



C 

CA 



CA 



n-l 



e A (t,0)[B AB 



A n ~ l B\ = 0. 



(5.4) 



But, the system is controllable and observable, and so we can form invertible ma- 
trices Sfp and &q by choosing n independent columns of the controllability matrix 
and n independent rows of the observability matrix, respectively. Then, by ([57 
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Hmt^oo^o e ^(*>0)^C = 0. Hence, \im t ^oo eA{t,0) — and exponential stability 
follows from the arguments given in Theorem 15.41 □ 



We use the preceding theorem in the following example. 
Example 5.14. Suppose T is a time scale with < /i < 4. The system 



- 8 


1 " 




"2" 


45 


3 °1 


x(t) + 


1 


. 45 


10. 





u(t), x(0) = x , 



y(t) = [3 i]x(t). 



is controllable by Example 12.81 and observable by Example 13.61 The eigenvalues of 
A are Ai = — § and A2 = — g. Note that the assumption on T implies Ai, A2 S 5(C), 
the stability region of T. Thus, by Theorem 1 5. 41 the system is exponentially stable. 
Theorem 15. 131 then says that the system is also BIBO stable. 

As we have seen, the Laplace transform can be a useful tool for analyzing stability 
in the time invariant case. With this in mind, we desire a theorem that determines if 
a system is BIBO stable by examining its transfer function. The following theorem 
accomplishes this. 

Theorem 5.15 (Transfer Function Characterization of BIBO Stability). The time 
invariant regressive linear system 

x A (t) = Ax(t) + Bu{t), x(t )=x , 

y(t) = Cx(t), 

is bounded-input, bounded- output stable if and only if all poles of the transfer func- 
tion G(z) = C{zl — A)~ 1 B are contained in S(C). 

Proof. If each entry of G{z) has poles that lie in 5(C), then the partial fraction 
decomposition of G(z) discussed earlier shows that each entry of G(t) is a sum of 
"polynomial-times-exponential" terms. Since the exponentials will all have sub- 
scripts in the stability region, 

||G(£)|| At < 00, (5.5) 





and so the system is bounded-input, bounded-output stable. 

Conversely, if (15. 5|) holds, then the exponential terms in any entry of G(t) must 
have subscripts in the stability region by using a standard contradiction argument. 
Thus, every entry of G(z) must have poles that lie in the stability region. □ 
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